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Abstract 

A Fokker-Planck type equation for interacting particles with exclusion principle 
is analysed. The nonlinear drift gives rise to mathematical difficulties in controlling 
moments of the distribution function. Assuming enough initial moments are finite, 
we can show the global existence of weak solutions for this problem. The natural 
associated entropy of the equation is the main tool to derive uniform in time a 
priori estimates for the kinetic energy and entropy. As a consequence, long-time 
asymptotics in are characterized by the Fermi-Dirac equilibrium with the same 
initial mass. This result is achieved without rate for any constructed global solution 
and with exponential rate due to entropy /entropy-dissipation arguments for initial 
data controlled by Fermi-Dirac distributions. Finally, initial data below radial solu- 
tions with suitable decay at infinity lead to solutions for which the relative entropy 
towards the Fermi-Dirac equilibrium is shown to converge to zero without decay 
rate. 

1 Introduction 

Kinetic equations for interacting particles with exclusion principle, such as fermions, have 
been introduced in the physics literature in P [121 [131 [111 [3 [23] and the review [TU] . 
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Spatially inhomogeneous equations appear from formal derivations of generalized Boltz- 
mann equations and Uehling-Uhlenbeck kinetic equations both for fermionic and bosonic 
particles. The most relevant questions related to these problems concern their long-time 
asymptotics and the rate of convergence towards global equilibrium if any. 

The spatially inhomogeneous situation has been recently studied in [22], where the 
long time asympotics of these models in the torus is shown to be given by spatially 
homogeneous equilibrium given by Fermi-Dirac distributions when the initial data is not 
far from equilibrium in a suitable Sobolev space. This nice result is based on techniques 
developed in previous works flUi fIT\ . Other related mathematical results for Boltzmann- 
type models have appeared in [H [19]. 

In this work, we focus on the global existence of solutions and the convergence of solu- 
tions towards global equilibrium in the spatially homogeneous case without any smallness 
assumption on the initial data. Preliminary results in the one-dimensional setting were 
reported in [3] . More precisely, we analyse in detail the following Fokker-Planck equation 
for fermions, see for instance [10], 



with initial condition f{0,v) = fo{v) G L^(M^), < /o < 1 and suitable moment condi- 
tions to be specified below. Here, f = f{t, v) is the density of particles with velocity v at 
time t > 0. 

This equation has been proposed in order to describe the dynamics of classical inter- 
acting particles, obeying the exclusion-inclusion principle in [12]. In fact, equation (11.11) 
is formally equivalent to 



from which it is easily seen that Fermi-Dirac distributions defined by 



with /3 > are stationary solutions. Moreover, for each value of M > 0, there exists a 
unique p = (3{M) > such that has mass M, that is, = M. Throughout 

the paper we shall denote by Fm- 

Another striking property of this equation is the existence of a formal Liapunov func- 
tional, related to the standard entropy functional for linear and nonlinear Fokker-Planck 
models [H [2] , given by 



dl 

at 



AJ + div,[vf{l-f)] 



V e w\t > 0, 



(1.1) 
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We will show that this functional plays the same role as the H-functional for the spatially 
homogeneous Boltzmann equation, see for instance [2lj. In particular it will be crucial to 
characterize long-time asymptotics of fll.ip . In fact, the entropy method will be the basis 
of the main results in this work; more precisely by taking the formal time derivative of 
we conclude that 



^W) = -/ /(I-/) 



V„ log 



f 



dv < 0. 



Therefore, to show the global equilibration of solutions to fll.ll) we need to find the right 
functional setting to show the entropy dissipation. Furthermore, if we succeed in relating 
functionally the entropy and the entropy dissipation, we will be able to give decay rates 
towards equilibrium. These are the main objectives of this work. Let us finally mention 
that these equations are of interest as typical examples of gradient flows with respect to 
euclidean Wasser stein distance of entropy functionals with nonlinear mobility, see [H |3] 
for other examples and related problems. 

In section 2, we will show the global existence of solutions for equation f ll.ip based on 
fixed point arguments, estimates involving moment bounds and the conservation of certain 
properties of the solutions. The suitable functional setting is reminiscent of the one used 
in equations sharing a similar structure and technical difficulties as those treated in [8], [11] . 
The main technical obstacle for the Fermi-Dirac-Fokker-Planck equation fll.ip lies in the 
control of moments. Next, in section 3, we show that the constructed solutions verify 
that the entropy is decreasing, and from that, we prove the convergence towards global 
equilibrium without rate. Again, here the uniform-in-time control of the second moment 
is crucial. Finally, we obtain an exponential rate of convergence towards equilibrium if 
the initial data is controlled by Fermi-Dirac distributions and the convergence to zero of 
the relative entropy when controlled by radial solutions. 



2 Global Existence of Solutions 

In this section, we will show the global existence of solutions to the Cauchy problem to 
fll.ip . We start by proving local existence of solutions together with a characterization 
of the time-span of these solutions. Later, we show further regularity properties of these 
solutions with the help of estimates on derivatives. Based on these estimates we can 
derive further properties of the solutions: conservation of mass, positivity, L°° bounds, 
comparison principle, moment estimates and entropy estimates. All of these uniform 
estimates allow us to show that solutions can be extended and thus exist for all times. 
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2.1 Local Existence 



We will prove the local existence and uniqueness of solutions using contraction-principle 
arguments as in [H El [H] for instance. As a first step, let us note that we can write fll.ll) 

as 



(2.1) 



and, due to Duhamel's formula, we are led to consider the corresponding integral equation 



f{t,v) 



J-'{t,v,w)fo{w)dw — / / ^{t — s,v,w){divw{wf{s,w) ))dwds (2.2) 

'o Jr'^ 



where JF(t, v,w) is the fundamental solution for the homogeneous Fokker-Planck equation: 

df 



dt 



div,{vf + VJ) 



given by 
with 



iV i 
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a{t) := e"2* , u{t) := e^* - 1 and Mi{^) := (27r/)~^e' 
for any A > 0. Let us define the operator J^{t,v)[g] acting on functions g as: 

J^{t,v)[g{w)]= I J='{t,v,w)g{w)dw. 



(2.3) 



Note that by integration by parts, the expression J^(t,v)[diVw{u!p{w))] is equivalent 



to: 



{271 (e2* - 1))^ 



{2n (e2* - 1))^ 
e~* (Vt,JF(t, ti, w) ■ w) f{wY dw 



f{wY dw 



so that 02.21) becomes 



f{t,v)=T{t,v)[U{w)]+ / e-(*-^)V.^(t-s,t;)K(s,«7)2]ds. (2.4) 

^0 
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We will now define a space in which the functional induced by (12.41) 



T[f]it,v):=J^it,v)[foiw)]+ / e~^'-''>V,Ht-s,v)[wf{s,wy]ds (2.5) 



has a fixed point. To this end, we define the spaces T := L°°(R^) n L}(M^) n Lp,(M^) 
and Tt '■= C([0,T]; T) with norms 

||/(t)||x := max{||/(t)|U, \\fit)\\^., \\fit)\\^.J and ||/||t, := max ||/(t)||T 

for any T > 0, where we omit the N-dimensional euclidean space for notational 
convenience and 



L.„:=||(l + |^;r)/||, and ||/||, := / l/^dt; 



In the following, we will see that for p > N, p > 2, and m > 1 we can choose q and r 
satisfying 

Np p mp P ^ ^ foR\ 
< - < r < < p and - < q < p 2.6 

such that !!'?'[/] IItt is bounded by ||/||tt- Let us fix such parameters p, m, r, q and 
< t < T. Due to Proposition lA.ll and q < p < 2q, we can compute 

pt Nit-s) 

||T[/](t)|U<Ce^*||/o||oo+ / C ^^|||^,|/2(,)||^d. 

rt pN{t-s) p p 

<Ce^*||/o||oo+ / C- -^||/(.)||oo"||/(s)||,%d. 

ft Nit-s) 

<Ce^*||/o||oo+ / C- ^ds 11/11^^ 

<C^e^*||/o||oo + CJi(t) 11/11^^, 



where 



by the choice (12.61) of q. In the same way, since r satisfies (m + l)r < mp and 2r > p, we 
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get 



<CeV'\\M\L?. + I C ^ .,,. \\fis)C"\\fis)\\ys 



where 



<Ce^'\\f,\\^.^+CUt) ll/lll^, 



by the choice (12.61) of r. 
Finally we can estimate 

.t (J 



\\mmLi<c\\foki+ I -—--rWHfisn^ids 

u[t — s)-i 





where by interpolation, we get as p > 2 and m > 1 

p2|| / /I , l„..IM„..l ^2 ^ / /I , \„..\\2i:2 



w\nLi= {i + \w\)\w\f'dw< {i + \w\yf'dw 

p-2 1 

<(/ {l + \w\)fdwY'{ f {l + \w\YfPdw 

<ii/iirii/ii;?- (2-7) 



Consequently 



\\mmLi<C\\fo\\Li+C I X-^(1-X)~^dx 11/11^, 



-2t 



We next check the existence of a fixed point of (12.51) in Tt- To this end, we define a 
sequence (/n)n.>i by = T[/„] for n > 0. Collecting all the above estimates, we can 
write 

||/n+i(t)||T<Ci(iV,t)||/o||T + C2(iV,p,g,r,t)||/„,f^^ 
for any < t < T and any T > 0, with 
Ci{N,t) := Ce^* 

C2(iV,p,g,r,t) := C max { I,{t),l2{t), [ X~H^ - x)~~'dx 

Je--2t 

6 



which are clearly increasing with t and C2{t) tends to as t does. Thus, for any T > 

||/„+i||tt <Ci(T) ||/o||^ + C2(T) ||/„|||^ 

with Ci(T) = Ci{N,T) and C2(T) = C2{N,p,q,r,T), both being increasing functions of 
T. We may also assume that Ci(T) > 1 without loss of generality. 

From now on, we will follow the arguments in [18]. We will first show that if T is small 
enough, the functional T is bounded in Ty, which will in turn imply the convergence. Let 
us take T > and 6 > which verify 

II/oIIt < S and < S < 



AC,{T)C2{Ty 

Then, let us prove by induction that ||/n||TT < 2Ci(T)(5 for all n. By the choice of T and 
5 we have ||/o||t < Ci{T)5 < 2Ci(T)S. If we suppose that ||/n||TT < 2Ci(T)5, we have 

ll/n+iIlT, < C^{T)5 + AC!{T)C2{T)5^ < 2C^{T)5, 

hence the claim. Now, computing the difference between two consecutive iterations of the 
functional and proceeding with the same estimates as above, we can see for any < t < T 
that 



\\fn+i — fi 



n Tt 



Jo 



< C2{T) sup \\fn + fn~l\\J\\fn " fn^l\\y^ 

[0,T] 

< C2{T) (^ll/nll^j, + ll/n-lllTr) H'^" ~ •^"-i|Itt 

< AC,{T)C2{T)5\\U - fn-i\\^^ < (4Ci(T)C2(T)5)"||/i - /o||^^ . 

Since 4Ci(T)C2(T)5 < 1 we can conclude that there exists a function in which is 
a fixed point for T, and hence a solution to the integral equation (12. 2p . It is not difficult 
to check that the solution / G to the integral equation is a solution of (11. ip in the 
sense of distributions defining our concept of solution. We summarize the results of this 
subsection in the following result. 

Theorem 2.1 (Local Existence) Let m>l, p>N,p>2, and /q G T. Then there 
exists T > depending only on the norm of the initial condition /o in T , such that (11.11) 
has a unique solution f in C([0,T]; T) with /(O) = /q. 

Remark 2.2 The previous theorem is also valid for /o G (L°° fl L^ fl L^)(]R^), with a 
solution defined in C([0,T]; (L°° nL^nL-'^)(]R^)) but we will need to have the first moment 
of the solution bounded in order to be able to extend it to a global in time solution. We 
thus include here this additional condition. 
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Remark 2.3 With the same arguments used to prove Theorem \2.1\ we can prove an equi- 
valent result for the Bose-Einstein-Fokker-Planck equation 

^ = AJ + dw^vfil + /)], veR'',t> 0. 
2.2 Estimates on Derivatives 

Let us now work on estimates on the derivatives. By taking the gradient in the integral 
equation, we obtain 

VJ{t,v) = V,J^{t,v)[f{w)]-f V,J^{t-s,v)[diY^{wf{s,w))]ds. (2.8) 

Jo 

where V vJ^{t,v)[g\ is defined as the vector: 

V^J^{t,v)[g\:= I V^J^{t,v,w)g{w)dw 

for the real-valued function g. Here, we will consider a space Xt with suitable weighted 
norms for the derivatives 

= max<^ ll/llxy, sup iy{t)^\\\/^f{t)\\-^P^, sup i^it)^\\^vfit)\y \ 

I 0<t<T 0<t<T ) 

where for notational simphcity we refer to || | V^/| Hl^^ as || Vt,/||Lf„. Let us estimate the L^- 
and L^-norms of V^,/ using again the results in Proposition lA.ll as follows: for r G 
satisfying (12. 6p 

i]t 



|V./(t)||LP^ < C^^-^||/o||l|;,+ / \\V,,J'[2f{wV^f)]+Nf\y^ds 

U{t)2 Jo 

< C*— -^||/o|ks. + C / — -^||/(s)||L^J|/(s)|Ud. 

iy{t)2 Jo u(t — s)'^ 

+C Z7r^ll/(^-V./)||L^ds 

Jo vit — s)^^^ p>^ 

Q\V' J n iV+2p' 1 



+C sup {K^)^/^||/(s)(^-V./(.))||Lj/(t) 

0<s<T 
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where 



-t /•! 



m < - 



1 f( N+r' 

- 2 



N + 2r 



1 + c 



1 -X) ^^^^ "^^^Hx-e 5dx 



1 - e- 



l+e- 



1 - e- 



-2t 



< Ce*^(l -e 

< Ch{t)v{t)-^2 



1 — K(1 — L\ 

2 2 W p' 



^/ N+p' ■, , l_iVM_].-, 

with := e p' V(t)2 2 "-r p^ which is an increasing function of time with /i(0) 
since p > r > Np/{N + p). It remains to estimate \\f{w ■ Vwf)\\L^- 



\w 



(m.+l)r rr 



Now, we can bound these integrals by using Holder's inequality to obtain 

p—r 

Jrn 

and 

\w 



(m+l)r rr 



nv^fi^dw < 



p — r 
P 



iwriv^f^dw 



Since p < pr /{p — r) < mp or equivalently (m + l)r/m < p < 2r by (12.61) . we have for 
any < t < T 



nV^fYdw < 



00 



riiv^/ii; < 



2r 



and 



k|(™+^)viv./r-d^< ii/iirii/rLriiv-/iiL^. < 



2r 



u{ty. 



Putting together the above estimates we have shown that, 

iy{tnm{wv^fm\Lr^<c\\f\\ 



9 



and 

^{t)H^vf{t)\U < Cl{T,N,p)\\fo\U+C',{T,N,p,r)\\f\\j,^ (2.9) 
with Cl and increasing functions of T and for any < t < T. Analogously, we reckon 

l|V./(t)|Li < C^II/oIIlj + C f /'"\j |/(.)||oo||/(^)||li ds 

+^ r ,f ^V ii/(^-v.,/)(.)iiLidg 

where by taking p > 2 and by interpolation as in (12.71) . we have 

||/(^-v./)||li < Ilkl^/||2|lkl^|v,./|||2 

< Il/ll|^ll/ll5^liv./||f^||v./||jf^ 

^ ll/llxr 
- ^^(t)V2' 

Putting together the last estimates, we deduce 

Ki)^l|V./(t)||Li <C3(T,iV,p)||/o||Li+C|(r,iV,p,r)||/||^^ (2.10) 

with Cf and C| increasing functions of T, for any < t < T. From (12.91) and (I2.10p and 
all the estimates of the previous section, we finally get 

ll/IU, <Ci(T,iV,p)||/o||T + C2(T,iV,p,r)||/||^^ 

for any T > 0. From these estimates and proceeding as at the end of the previous section, 
it is easy to show that we have uniform estimates in Xt of the iteration sequence and the 
convergence of the iteration sequence in the space Xx- From the uniqueness obtained in 
the previous section, we conclude that the solution obtained in this new procedure is the 
same as before and lies in Xj-- Summarizing, we have shown: 

Theorem 2.4 Let m>l,p>N,p>2, and fo G T. Then there exists T > depending 
only on the norm of the initial condition /o in T such that (II. ip has a unique solution 
in C([0,T];T) with /(O) = fo and velocity gradients verifying that t ^ i^(t) 2 1 V„/(t)| G 
fiC((0,T),(L^nLi)(M^)). 

2.3 Properties of the solutions 

As (11.11) belongs to the general class of convection-diffusion equation, it enjoys several 
classical properties which we gather in this section. The proof of these results uses classical 
approximation arguments, see [HI [2S] for instance. Since these arguments are somehow 
standard we will only give the detailed proof of the L^-contraction property below. 
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Lemma 2.5 (Positivity and Boundedness) Let f G Xt be the solution of the Cauchy 
problem fll.ip with initial condition /o G T. If < Jq < 1 in M^, then < f{t) < 1 for 
any < t < T . 

Lemma 2.6 (L^-Contraction and Comparison Principle) Let f G Xt and g G Xt 

be the solutions of the Cauchy problem (11.11) with respective initial data /o G T and g^ G T. 
Then 

||/(t)-^?(t)||i< ||/o-^7o||i (2.11) 
for all < t < T . Furthermore, if /o < go then f{t,v) < g{t,v) for all < t < T and 
V G M^. 

Proof.- Since / and g solve fll.ll) . 

l(f - g) = A^if - g) + VMf - 9)) - ^Mf - /)) (2.12) 

holds. We will obtain this result from the time evolution of | / — g\i; where | ■ |e is the 
primitive vanishing at zero of sign£(s), the latter being an increasing smooth approxima- 
tion of the sign function defined by sign (s) = 1 if s > 0, sign (0) = and sign (s) = — 1 
if s < 0. Multiplying both sides of equation (12.121) by Cn(^) sign£(/ — g) and integrating 
over R^, where £ C^(]R^) is a cut-off function satisfying < < 1, Cn(^) = 1 if 
1^1 < ''^5 Cn(^) = if |f I > 2n, and |V„Cn| < ^5 we obtain 

^/ U^)\f-9\edv<-[ Uv)signi{f-g){vVM-9)){f-9)dv 

+ [ Uv) signU/ -g){v- V.(/ - g)){f - /) dv 

- / V.CnSign,{f-g){VM-9) + vif-9-if-9')))dv 

= - [ Cn{v){v ■ V.((/ - g) sign,(/ - g) - \f - gU)) dv 
Jrn 

+ I Uv){f + 9){v - W^iif - g)sign,{f - g) - \f - gl))dv 
Jrn 

- [ V,CnSign,{f-g)iV.„if-g) + v{f-g-{f-g')))dv. 

JRN 



Integrating by parts, we finally get 
dt 



4/ Cn{,v)\f - g\edv < I div^«„,(i;))((/-^)signe(/-^) - |/-^ydw 



- f diY,,{Uv)v{f + gMf - g) sign, {f - g) -\f - gl)dv 
Jrn 

+ -/ \'^v{f-9) + vif-9-{f-9'))\dv. 
n Jrn 
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For every n, the first two integrals become zero as e — > 0, since / and g are in Xt whence 
fit), git) G n L°°(M^) and V„/(t), V„^(t) G Lj(M^) for any < t < T, allowing for 
a Lebesgue dominated convergence argument. We have that V^,/ + vf{l — /) G L^(R^) 
and V^g + ^^'(l — g) E L^(]R^) for any < t < T, and thus the third integral vanishes as 
n — > oo, getting finally 

\f-g\dv<0 (2.13) 

which concludes the proof of the first assertion of the lemma. □ 
Similar arguments show the conservation of mass. 

Lemma 2.7 (Mass Conservation) Let f G Xt be the solution of the Cauchy problem 
(11. ip with non-negative initial condition fo G T, then the L^-norm of f is conserved, i.e. 
ll/Wlli = ll/o||i/or allte[0,T]. 

Finally, we establish time dependent bounds on moments of the solution to (II. ip . More 
precisely, we will show that moments increase at most as a polynomial on t. First, let us 
note that given a, 6 > 1 and / G L^ft(M^) n L°°(M^) then 

II/IIU<C||/||^ ||/||L"^ (2.14) 

ab 

Indeed, 

ii/iiL.=f/ {i+ivirfdvY <(c f (i+i^nfd^y 

< (cwfif / (1 + \vnfdv) " = c\\f wi wftJ. 

In particular, (L^^ n L°°)(M^) C T. 

We next define [7] to be the smallest integer larger or equal than 7. 

Lemma 2.8 (Moments Bound) Let f G Xt be the solution of the Cauchy problem 
(II. ip with initial condition /o G L^^p(M^) for some m>l,p>N,p>2, and satisfying 
< /o < 1- Then, for <t <T and 1 < 7 < mp/2 the 2'~f-moment of fit) is bounded by 
a polynomial P^^-T^it) of degree [7], which depends only on the moments of f^. 

Proof.- We will prove it by induction on 7. First, we will see that the second moment is 
bounded, and afterward that we can bound every moment of order smaller than pm in 
terms of a 7*^^ moment with < 7* < 2, which can in turn be bounded in terms of the 
second moment. 
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Let (Cn)n>i be a sequence of smooth cut-off functions satisfying < Cn. < 1, Cn{v) = 1 
if l^^l < n, Cn{v) = if \v\ > 2n, |V^Cn| < and |A^Cn| < V^^- We multiply ([HI]) by 
I^PCn(^) and integrate over to get 

d_ 
dt 



Uv)\v\'fdv= Uv)\v\'AJdv+ Uv)\v\Mw,{vf{l-f))dv 

< f [A,Cn\v\^ + 4V.CnV + 2NCn]fdv+ [ | V.Cn | | '/( 1 " /)dw 



2 / Cn\v\'fdv + 2 / CnK^rrd^ 



<5 / fdv + 2N Cnfdv+ / l^r/dw. 

-'n<|t)|<2n JM^ >/n<|«;|<2n 

Now, letting n ^ oo and noticing that fl{n<\v\<2n} and |f P/l{n.<|D|<2n.} converge pointwise 
to zero and are bounded by / and \v\^f respectively with / G X^, we infer from the 
Lebesgue dominated convergence theorem that the first and the last integrals converge to 
zero. Finally, integrating in time, we get 

\v\^f{t,v)dv<l \v\^ fQ{v)dv + 2N Mt (2.15) 



for all < t < T. Now, for the moment 27 we can see in the same way 
di 



Uv)\v\'ydv= I Uv)\v\'^AJ-dv+ I Uv)\v\'Miv,{vf{l-f))dv 



< / [A,Cn\v\'^ + A-fV,Cn\v\'^^-'^V + 27(2(7 - 1) + N)\v\^^^-'^Cn] fdv 

Jrn 

+ [ |V,C„||^|'"+V(l-/)dt'-27 / Cn|^|'Vdt' + 27 / Cn\v\''fdv 



<C |tf (^^^Vdf + 27(2(7 - 1) + AT) / Cn|fr^^"'Vdw 

Jn<\v\<2n JrI^ 
J n<\v\<2n 

and we again let n go to infinity. If 27 < mp, the previous argument ensures that only 
the second integral remains, and integrating in time, we conclude 

\v\^y{t,v)dv < [ |t;|2^/o(t;)d^; + 27(2(7 - 1) + iV) / / \v\^^^-^^ f{s,v)dvds 

(2.16) 

for all < t < T. Whence, if we assume by induction that the hypothesis of the lemma 
holds true for the 2(7 — l)-moment, 

\v\'y{v,t)dv< [ \v\'yo{v)dv + 2^{2{^-l) + N) fp^,_^^{s)ds (2.17) 
Jr^ Jo 
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for all < t < T, defining by induction the polynomial P|-^-] . □ 
2.4 Global existence 

Given an initial condition /o G L^p(]R^), p>N,p>2, m>l such that < /o < 1, 
we have /o G T and we have shown in the previous subsections that there exists a unique 
local solution of (11. ip on an interval [0,T). In fact, we can extend this solution to be 
global in time. If there exists T^ax < oo such that the solution does not exist out of 
(0, Tmax), then the T-norm of it shall go to infinity as t goes to T^ax] as "we will see, that 
situation cannot happen. 

Due to Lemma [2.51 we have that < f{t, v) < 1 for any < t < T and any v G M^, 
and thus a bound for the L°°-norm of f{t). Also, the conservation of the mass in Lemma 
12.71 together with the positivity in Lemma 12.51 provide us with a bound for the L^-norm. 
Finally, due to (12.141) and Lemma 12.81 the L^-norm is also bounded on any finite time 
interval. 

Theorem 2.9 (Global Existence) Let fo G L;^p(R^), p>N,p>2, m>lbe such 
that < /o < 1. Then the Cauchy problem (11.11) with initial condition fo has a unique 
solution defined in [0, oo) belonging to Xj- for all T > 0. ^4/50, we have < f{t,v) < 1, 
for allt>0 andv G and ||/(t)||i = ||/o||i = M for all t > 0. 

Remark 2.10 Note that for any K > we can consider (11.11) restricted to the cylinder 
Ck '■= [0, oo) X {|t;| < K}. Then, due to the fact that the solutions to (11.11) we have con- 
structed are in L°° , we can show that the solution is indeed C°°{Ck) by applying regularity 
results in [16] for quasilinear parabolic equations. 



Corollary 2.11 If fo G L^p(]R^) fl L°°(M^) is a radially symmetric and non-increasing 
function (that is, fo{v) = ipo{\v\) for some non-increasing function ipo), then so is f{t) 
for all t > 0, that is, f{t, v) = (p{t, \v\) and r i— > if(t^ r) is non-increasing for all t > 0. In 
addition, ip solves 

' ^ ^r^-^+rVl-^)! n^^th ^{t,0)=0 (2.18) 



dt r^~i dr \ dr J dr 

and ip{0, r) = foi^)- 

Proof.- The uniqueness part of Theorem 12.91 and the rotational invariance of (II. ip imply 
that f{t) is radially symmetric for all t > 0. The other properties are proved by classi- 
cal arguments, the monotonicity of r i— > (p(t, r) being a consequence of the comparison 
principle applied to the equation solved hj d(p/dr. □ 
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3 Asymptotic Behaviour 



Now that we have shown that under the appropriate assumptions equation f 1 1.11) have a 
unique solution which is global in time, we are interested in how does this solution behave 
when the time is large. For that we will define an appropriate entropy functional for the 
solution and study its properties. 



3.1 Associated Entropy Functional 

In this section, we will show that the solutions constructed above satisfy an additional 
dissipation property, the entropy decay. For (? G T such that < g < 1, we define the 
functional 

Hig):=S{g) + Eig) (3.1) 

with the entropy given by 

S{g):= [ s{g{v))dv (3.2) 



where 

s(r) := (1 -r)log(l -r) +rlog(r) < 0, re[0,l], (3.3) 
and the kinetic energy given by 

E{g):=l- [ \v\'g{v)dv. (3.4) 
^ Jrn 

We first check that H{g) is indeed well defined and establish a control of the entropy in 
terms of the kinetic energy. 

Lemma 3.1 (Entropy Control) For e G (0, 1), there exists a positive constant such 
that 

< -S{g) < eE{g) + C, (3.5) 
for every g G L2(M^) such that < g < 1. 

Proof.- For e G (0, 1) and v G M^, we put Zs{v) := 1/(1 + e^''''^^^). The convexity of s 
ensures that 

- s{ze{v)) > s'{z,{v)){g{v) - Ze{v)) 
-s{z,{v)) + s{g{v)) > log (y^II^) iaiv) - z,{v)) 
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for V e M^. Since z^{v)/{l — Ze{v)) = e ^I'^l^/^, we end up with 
-s{g{v)) < -^g{v) - s{z,{v)) —Ze{v) 



2 

e\v\ 



g{v) + (1 - z,{v)) log (l + e-^l^'l'/^) + bg (l + g-^l'^l'/^) 



2 

< + e-HV2 (3.6) 

for G M^, where we used log(l + a) < a for a > and < < 1. Integrating the 
previous inequahty yields (13. 5p . □ 

We next recall that Fm is the unique Fermi-Dirac equilibrium state satisfying H-Fa/Hi = 
M := II /oil i; then we can introduce the next property for H. 

Lemma 3.2 (Entropy Monotonicity) Assume that f is the solution to the Cauchy 
problem (II. ip with initial condition /o in L^p(M^) for some p > max{N,2) , m > 1 and 
satisfying < /o < 1 . Then, the function H is a non-increasing function of time satisfying 
for allt>0 that 

H{fo)>H{f{t))>H{FM) with M:=||/o||i. (3.7) 

Proof. - We first give a formal proof of the time monotonicity of H{f) and supply additional 
details at the end of the proof. First of all, we observe that we can formulate (II. ip as 



A- 



/(I - /)V. ( s'if) + ^ 



2 

We multiply the previous equation by s'{f) -\-\v\'^/2 and integrate over to obtain that 

^H{f) = -f f{l-f)\v + Vj{f)\'dv<Q. (3.8) 
at J^N 

Consequently, the function t — > H{f(t)) is a non- increasing function of time, whence 
the first inequality in (13. 7p . To prove the second inequality, we observe that the convexity 
of s entails that 

s{f{t,v)) - s{Fm{v)) > s'{FM{v)){f{t,v) - Fm{v)) 

s{Fm{v)) - s{f{t,v)) < (\og(3{M) + {f{t,v) - Fm{v)) 

for {t,v) G [0, oo) X M^. The second inequality in (13.70 now follows from the integration 
of the previous inequality over since ||-Fa/||i = ||/(^)||i by Lemma [27fl 
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We shall point out that, in order to justify the previous computations leading to 
the time monotonicity of the entropy, one should first start with an initial condition /, 
e G (0, 1), given by 



0; 



max < mm 



/o(tO, 



1 



eel 



V G 



Owing to the comparison principle (Lemma 12.61) . the corresponding solution to (11. II) 
satisfies 



< ^ < nt,v) < ^ , \ ,,,, < 1, {t,v) G (0,oo) X 



e + e 



1 +5el-l'/2 



for which the previous computations can be performed since the solutions are immediately 
smooth and fast decaying at infinity for all t > 0, and thus if(/^(t)) < H^Jq) for all t > 0. 



Since /q — * /q in T and in 



mp \ 



as e 0, it is not difficult to see that redoing 
all estimates in subsections 2.1 and 2.2, we have continuous dependence of solutions with 
respect to the initial data, and thus, converges towards / in Xt for any T > 0. 
Moreover, we have uniform bounds with respect to e of the moments in finite time intervals 
using Lemma YIM Direct estimates easily show that H{f^) — > H^Jq) as e — > 0. 

Let us now prove that H{f^(t)) H{f(t)) as e — > for t > 0. Let us fix i? > 0. Since 
/^(t) f{t) in L^(]R^) and we have uniform estimates in e of moments of order mp > 2 
then 



\v\\nt)-f{t))dv 



< 



< 



\v\'\nt)-f{t)\dv 



v\>R 
1 



\v\\nt)-f{t))dv 



'\v\<R 

r{nt) + fit))dv 



v\>R 



< 



+R'\\nt)-fit)\\, 

C{t) 



R'\\nt)-m\\, 



Since the above inequality is valid for all i? > 0, we conclude that E{f^{t)) E{f(t)) 
as e -> 0. Now, taking into account that (1 + \v\'^)Pit) (1 + \v\'^)fit) in L^(M^), we 
deduce that there exists h G L^(]R^) such that ||fp/^(t)| < h and f^(t) f{t) a.e. in 
M^, for a subsequence that we denote with the same index. Using inequality (13. 6p . we 
deduce that 

< ~s{f{t,v)) < ^/i(i;) + e-l"l'/^ G LH 



and that —s{f'^(t,v)) — s(/(t,t>)) a.e. in M^. Thus, by the Lebesgue dominated 
convergence theorem, we finally deduce that S{f^{t)) S{f(t)) as e — >• 0. The con- 
vergence as e — > of S'(/^(t)) to S{f{t)) is actually true for the whole family (and not 
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only for a subsequence) thanks to the uniqueness of the hmit. As a consequence, we 
showed H^f^it)) H{f(t)) as e ^ and passing to the hmit e ^ in the inequahty 
Hifit)) < H{f^), we get the desired result. □ 

Now, it is easy to see the existence of a uniform in time bound for the kinetic energy 
E{f(t)), or equivalently, of the solutions in L2(M^). If we take equations (13.11) . (13.51) (with 
e = 1/2) and (13.71) we get that 

E{f{t)) = H{f{t)) - S{f{t)) < lE{f{t)) + Cy, + H{fo) 
for t > whence 

E{fit))<2{C,/2 + Hifo)). (3.9) 



3.2 Convergence to the Steady State 

Theorem 3.3 (Convergence) Let f be the solution to the Cauchy problem (II. ip with 
initial condition fo in L^p(]R^), p > max{N,2), m > 1 satisfying < /o < 1. Then 
{f(t)}t>Q converges strongly in L^(R^) towards Fm as t ^ oo with M := ||/o||i- 



For the proof, we first need a technical lemma. 



Lemma 3.4 Let f be the solution to the Cauchy problem (11.11) with initial condition Jq 
in L^p(M^), p > max{N,2), m>l satisfying < /o < 1. If A is a measurable subset of 
M^, we have 



I (I \vf{l - /) + V,/| dv\ dt < H{Fm) sup ( / f{t,v)dv\ 
Jo \Ja / t>o Ua ) 



(3.10) 



Proof.- Owing to the second inequality in (13. 7p and the finiteness of iJ(/o), we also infer 
from dSIHD that {t,v) i — > /(I - /) |v + V„s'(/)|^ belongs to L^{{0,oo) x M^). Working 
again with the regularized solutions it then follows from Lemma [2. 71 and the Cauchy- 
Schwarz inequality that, if A is a measurable subset of M^, we can compute 

- rf/ '-'^^^^-^^^+y^'(/-(i -/-))-/- d^Vdt 

Jo \Ja (/^(l - f^)f' J 

-dv) (^Jj'il-f')dv^ dt, 



\vni-fn + v,r 
/^(i - /^) 
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and thus. 



<sup( / rit,v)dv] r [ f{i-f)[v+w,s'{nfdvdt 

t>o Ua J Jo J a 

<i/(FMOsup( / r{t,v)dv\ . 

t>o Ua ) 

Here, := ||/o||i so that Fm^ is the Fermi-Dirac distribution with the mass of the 
regularized initial condition /q. It is easy to check that H{Fm^) H{Fm) as e ^ since 
—>■ M as e ^ 0. Passing to the limit as e ^ 0, ^ / in for any T > 0, and thus 
we get the conclusion. □ 

Proof of Theorem \3.3[ - We first establish that 

{f{t)}t>o is bounded in L^(M^) n L~(M^) . (3.11) 

From (13. 9p and Theorem 12. 9[ it is straightforward that E{f{t)) is bounded in [0,oo). 
Recalling the mass conservation, the boundedness of {/(t)}t>o in L2(M^)nL°°(]R^) follows. 

We next turn to the strong compactness of {f{t)}t>o in L^(M^). For that purpose, we 
put R(t,v) := vf(t,v){l — f{t,v)) for (t,v) G (0, oo) x and deduce from Theorem 12.91 
and dXTT]) that 

sup(||i?(t)||i + ||i?(t)||^) <2 sup /" {l + \v\^)f{t,v)dv <oo. (3.12) 

Denoting the linear heat semigroup on by (e*^)t>o, it follows from ([H]) that / IS given 
by the Duhamel formula 

/o+ / V^e(*-")^i?(s)ds, t>0. (3.13) 
Jo 



fit) = e 



It is straightforward to check by direct Fourier transform techniques that 

\\e'^9\\H^<C{a) min{r"/l^?|h,t-(^"+^)/1|^?||i} 
for teiO,oo), g e L\R^) n L2(M^) and a E [0, 2] with 

1/2 



\9\\h^--=([ ieP"i?(oi'rfA 
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and ^ being the Fourier transform of g. Thus, we deduce from (13.131) that, if t > 1 and 
a G ((1 - (A^/2))+,l), we have 



||/(t)b. < C(a)r(2-+^)/^||/o||i + C(a+l) / (t-.)-(2+2a+7V)/4||^^^)||^^^ 





+ C{a + l) [ (t-s)-(i+")/2p(s)||2ds 

< c, 

thanks to the choice of a. Consequently, {f(t)}t>i is also bounded in for a G 
((1 — (A^/2))"'", 1). Owing to the compactness of the embedding of (ij° fl L2)(M^) in 
L^(]R^), we finally conclude that 

{f{t)}t>o is relatively compact in L^(M^) . (3.14) 

Consider now a sequence {tn}nGN of positive real numbers such that t„ — *■ cxd as n — > oo. 
Owing to (13.141) . there are a subsequence of (not relabelled) and goo G L^(]R^) such 
that {f(tn)}nm converges towards goo in L^(M^) as n ^ oo. Putting fn(t) = f(tn + t), 
t G [0, 1] and denoting by g the unique solution to (II. ip with initial datum ^foo, we infer 
from the contraction property (12.111) that 

lim sup ||/„,(t)-(?(t)||i = 0. (3.15) 
"-^°°tg[o,i] 

Next, on one hand, we deduce from the proof of Lemma [3.41 with A = that {t, v) i — > 
vf{t, v){l - f{t, v)) + VJ{t, v) belongs to L\{0, oo); Li(M^)). Since 

t ( [ \vfn{l - fn) + dt = ( I \vf{l - /) + V,/| dv\ dt , 

Jo \Jrn J Jt^ \Jm J 

we end up with 

lim / (l |t;/„,(l-/n) + V,/„|d^;) dt = 0. (3.16) 

"^^•^ Jo \Jm / 

On the other hand, it follows from the mass conservation and (13.101) that, if A is a 
measurable subset of with finite measure \A\, we have 

1 / . X2 

vfn{l-fn)+Vvfn\dv] dt < if (Fm) | A | , 



\JA 
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which imphes that {vfn{l — fn) + ^vfn}nm is weakly relatively compact in L^((0, 1) x M^) 
by the Dunford-Pettis theorem. Since {t'/n(l — /n)}neN converges strongly towards vg{l — 
g) in L^((0, 1) x M^) by (13. Ill) and fl3.15p . we conclude that {V^/„}„>o is weakly relatively 
compact in -^^^((0, 1) x M^). Upon extracting a further subsequence, we may thus assume 
that {V^/n}n>o converges weakly towards VvQ in L-'^((0, 1) x M^). Consequently, 

/ / \vg{l- g) + S/,g\^vd.t<ymviTii f [ |t;/„(l - /J + dw dt = 

Jo Jm Jo 

by (13.161) . from which we readily deduce that vg{l — g) + V^g = a.e. in (0, 1) x . 
Since = M for each t G [0, 1] by Lemma I^TTI and (I3.15p . standard arguments allow 

us to conclude that g{t) = Fm for each t G [0, 1]. We have thus proved that Fm is the 
only possible cluster point in L^(M^) of {f(t)}t>o as t — >• oo, which, together with the 
relative compactness of {f(t)}t>o in L^(M^), implies the assertion of Theorem 13.31 □ 
By now, we have seen that the solution of (II. ip with initial condition /o converges to 
the Fermi-Dirac distribution Fm with the same mass as /o as t — > cx), but we are also 
interested in how fast this happens. We will answer that question with the next result, 
which was already proved in [5] in the one dimensional case, and easily extends to any 
dimension based on the existence and entropy decay results established above. 

Theorem 3.5 (Entropy Decay Rate) Let f be the solution to the Cauchy problem 
(II. ip with initial condition fo in L^p(]R^), p > max(iV, 2), m > 1 satisfying < /o < 
Fm* < 1 for some M* . Then 

H{f{t)) - H{Fm) < (Hifo) - H{FM))e-''^' (3.17) 

and 

||/(t) - FmWi < C^mh) - HiFM))'/'e-^' (3.18) 
for all t > 0, where C depends on M* and M := ||/o||i- 

Proof.- Since < /o < Fm*, then the initial condition satisfies all the hypotheses of 
Theorems 12.91 and 13.31 In order to show the exponential convergence, we use the same 
arguments as in [5]. We first remark that the entropy functional H coincides with the one 
introduced in [2] for the nonlinear diffusion equation 

^ = div^ [gV Ax + h{g))] (3.19) 

for the function < g(t,x) < 1, x G M, t > 0, where h{g) = s'{g) = logg — log(l — g). 
Let us point out that the relation between the entropy dissipation for the solutions of the 
nonlinear diffusion equation (I3.19p . given by 



-^o(^) = ^H{g) = - I 



r) ^ 
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and the entropy dissipation for the solutions of (11.11) . given by (13.81) . is the basic idea of 
the proof. Indeed, one can check that, once restricted to the range / G (0, 1), h{f) verifies 
the hypotheses of the Generahzed Logarithmic Sobolev Inequahty [21 Theorem 17]. The 
Generahzed Logarithmic Sobolev Inequality then asserts that 

Hig) - H{Fm) < \D,{g) (3.20) 

for all integrable positive g with mass M for which the right-hand side is well-defined 
and finite. We can now, by the same regularization argument as before, compare the 
entropy dissipation D{f) = — ^if(/) of equation (11.11) and the one -Do(/) of equation 
flCTj) . Thanks to LemmalSIHlwe have f{t,v) < Fm*{v) < {P{M*) + 1)"^ a.e. in M^, and 
thus 

D{f)= I f{l-f)\v + \/Mf)fdv>c[ f\v + VMf)fdv (3.21) 

where C = 1 — (/5(M*) + 1)~^. Applying the Generalized Logarithmic Sobolev Inequality 
(I3.20p to the solution / and taking into account the previous estimates, we conclude 

H{f{t)) - H{Fm) < {2Cr'D{f{t)). (3.22) 

Finally, coming back to the entropy evolution: 

I [H{f{t)) - H{Fm)] = -Difit)) < -2C [Him) - H{Fm)] , 

and the result follows from Gronwall's lemma. The convergence in L^ is obtained by a 
Csiszar-Kullback type inequality proven in [5], Corollary 4.3], its proof being valid for any 
space dimension. It is actually a consequence of a direct application of the Taylor theorem 
to the relative entropy H{f) — H{Fm) giving: 

\\f-FM\\l<2M{H{f)-H{FM)). 

□ 

3.3 Propagation of Moments and Consequences 

There is a large gap between Theorem 13.31 which only provides the L^-convergence to the 
equilibrium and Theorem 13.51 which warrants an exponential decay to zero of the relative 
entropy for a restrictive class of initial data. This last section is devoted to an intermediate 
result where we prove the convergence to zero of the relative entropy but without a rate 
for a larger class of initial data than in Theorem 13. 5[ 
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Lemma 3.6 (Time independent bound for Moments) Let go E L^p(]R^) withm > 
1, p > max [p, 2) such that < Qq < 1, and assume further that Qq is a radially symmetric 
and non-increasing function, i.e., there is a non-increasing function ipo such that go{v) = 
ipQ{\v\) for V G M^. Then, for the unique solution g of the Cauchy problem fll.ip with 
initial condition go, the control of moments propagates in time, i.e., there exists C > 
depending on N and go, but not on time, such that 

lim sup / \vrPg{t,v)dv = 0. (3.23) 

i>0 J{\v\>R} 

Proof.- We have already seen in Corollary 12 . 1 1 1 the existence and uniqueness of g and that 
g{t,v) = ip{t, \v\) for t > and v G for some function ip such that r i— (p{t,r) is 
non-increasing. Furthermore, we have that its moments are given by 

M := [ g{t,v)dv = Nlun [ r^-^ip{t,r) dr (3.24) 
Jm.^ Jo 

and 

\vr''g{t, v) dv = Nuon / r^+^P- V(t, r) dr (3.25) 

Jo 

for t > 0, where uj^ denotes the volume of the unit ball of M^. 

Next, since \v\'''^Pgo E L^(]R^), the map v h-* \v\"^p belongs to (R^ ; go{v) dv) and 
a refined version of de la Vallee-Poussin theorem P, [TTj ensures that there is a non- 
decreasing, non-negative and convex function i/j E C°°{[0,oo)) such that '0(0) = 0, ip' is 
concave, 

lim^^ = oo and / ip{\vrP)go{v) dv < oo. (3.26) 

Observe that, since ip{0) = and ip'{0) > 0, the convexity of tp and the concavity of ijj' 
ensure that for r > 

rilj"{r) < ifj'ir) and ipir) < rip'{r). (3.27) 
Then, after integration by parts, it follows from (12.181) that 

— i- [ ^(r'^Hr^-V dr = - [ r*"^" V(r"*^') ( r^-^^ + r^mh - m)] dr 
mpdtJo Jo \ dr J 

= h + h, (3.28) 

where 

/•oo 

h = [(mp + iV-2)r"^^+^-V(r"*^') +mpr2™^+^-V(r"*^')] dr 

Jo 

POD 

/2 = - / r^+™P" V(^"^)<^(1 - ^) dr. 
Jo 
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We now fix i? > such that unR^ > 4M and R'^ > A{2mp + iV - 2), and note that 
due to the monotonicity of ip with respect to r and fl3.24p - fl3.25p the inequality 

M > Nujm / r^~Vdr > ujnR^v{R) (3.29) 
Jo 

holds. Therefore, we first use the monotonicity of ip' and ip together with fl3.29l) to obtain 



h < - r^+"'P-^ij'{r'^PMl-ip)dr<{ip{R)-l) r^+^f" Vlr'^P)'^ dr 
Jr Jr 

V^A^-R J Jr ^ Jr 

o pR q poo 



4iVu;Ar 4 Jo 

On the other hand, from ((321, (ESS]), (E2ZD, (ESI) and the monotonicity of f 



oo 



h < {N + 2m-2) r^+'"^'-V(r'"^)(/9dr 



< (iV + 2mp - 2)^'(/2"P)/2'"P~2 /" 

io 



r (/? dr 



i?2 



oo 



Nun 4 ' -/^ v • 



Inserting these bounds for Ji and I2 in f l3.28p and using (13.271) we end up with 

1 rl r°° 

/ ^(r"^P)r^-Vdr 
mp at Jq 

< ^'(^""2^^"^''' (^ + N + 2mp-2)-- r r^+-f- V'(r"^^)<^ dr 
Nujn \ 4 J 2 Jq 



00 



iYcuAT \ 4 J 2 Jq 

We then use the Gronwall lemma to conclude that there exists C > depending on 
A^, M, m, p, and ijj such that 

sup / ^{\v\^'P)g{t,v)dv < C 



from which fl3.23p readily follows by fl3.26p . 
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Theorem 3.7 (Entropy Convergence) Let f be the solution of the Cauchy problem 
(11. ip with initial condition /o G L^p(]R^) such that there exists a radially symmetric and 
non-increasing function go G L^p(R^) with < fo < go < 1. Then H{f) H{Fm) as 
t ^ oo where M = ||/o||i- 

Proof.- Due to [191 Theorem 3] we know that 

\Hifit))~HiFM)\<C [ \v\'\f{t,v)-Fiv)\dw 

Jrn 

<R^f{t)-F\\, + snp [ \v\'\f{t)-F\dv 

t>0 J\v\>R 

Now, Theorem 13.51 and Lemma [3.61 imply that H{f(t)) H{Fm) as t — >■ oo. □ 



A LJ^-bounds for the Fokker-Planck Operator 

Here we follow similar arguments as in [H] to show some bounds for ||c^a^/(t)||L^ which 
were useful in the fixed point argument in Section 12.11 We recall the well-known Young 
inequality: Let gi G L'~(]R^), g2 G L'?(]R^) with 1 < p, r, g < oo and i + 1 = i + i, then 

5'i * 5'2 e L^(IR^) and US'! * 5'2||p < ||5'i||r ||5'2||g- 
Proposition A.l Let 1 < q < p < oo, m > and a G N^. Then for t > 0, 

C 



Proof.- For all a G N^, we have 



d^J'{t,v)[f] =d- ^e"^(^ f{w) dw 

Jr^ V(27r(e2* - 1))^ 



e 



d" / ^e f{e'w) dw 

Jr^ V(27r(e2* - 1))- 



gt(27V+|a|) , I y 



where 



f{e'w) dw (A.2) 



<l>M = d^i(l>o) (x) = ^h(x)0o(x), 
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being V\a\{x) ^ polynomial of degree |a| which we can recursively reckon by 



Voix) = 1, V\a\{x) = ^iahi(x) - X^|ahi(x) and Mx) = {27i)"e-—. 
Since 1 + < C{1 + \v - w|™)(l + we deduce 



(l+|i;r)ia^*/)(t)| < 



< C- 



^t{2N+\a\) 



N+\a 



p{t)^2— JrN 



(1 + Iv-wD 



V — w \ 



Then, we can write 



(1 + \v — w\ 



00 



V — w \ 
e-V(t)i/2 ) 



{l + lwD f{e'w) dw. (A.3) 



dw = C{I + //) 



with 



V — W \ I V — w 

00 



°l ^ e-V(t)V2 J 



e-V(t)i/5 



and 



II 



dw; 



mr)/2 



V — w \ , f V — w Y 



\xr'pl^iix)Mxr = C2 



whence 



On the other hand, we get 



(1 + \v - wry 



V — w \ 



e ''u 



dw < C. 



(A.4) 



fie'w) 



(1 + l^rf fie'w) 



dw 



e-^*(l + |e-*xrr /(X) 



dw 



dw 



(A.5) 



Putting (1A.4P together with (lA.Sp . we can use Young's inequality in (1A.3P as before. 



since 1 < q < p with r given by^ + l = ^ + ^to get the desired bound. 
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